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A MICROMETEOROID VELOCITY DETECTOR* 

By Frank Neuman 

SUMMARY 

A feasibility study for a particular type of particle velocity 
detector is made. The detector considered consists of two light 
screens with known spacing. A particle penetrating the successive 
light screens scatters light which is detected by means of photo- 
multiplier tubes . The time difference between the photomultiplier 
tube outputs is a measure of the particle velocity. 

Such detectors have been used many times in free-flight wind 
tunnels. The question that is answered here is whether such a 
system is useful In space for detecting extremely small particles 
which have very high velocities and which arrive from unknown direc- 
tions at random times with an unknown average rate. Can meaningful 
measurements be made in the presence of noise due to stray light? 

The answer is a 'qualified "yes." 


*A thesis submitted in June 1963 to the Department of Electrical 
Engineering and the Committee on the Graduate Division of Stanford 
University in partial fulfillment of the requirements for the degree 
of Engineer. 





INTRODUCTION 


The hazard which micrometeoroids present to space vehicles was 
recognized long before the first man-made satellites were flown. 

For this reason all the early American and Russian satellites carried 
some type of meteoroid detectors (ref. 1). At present, with manned 
space flight in mind, a satellite has been planned for the sole pur- 
pose of obtaining more information on the particle environment in 
space (ref. 2). One can think of two types of experiments. In the 
first type, which may be called engineering experiments, one exposes 
surfaces of important structural materials such as aluminum and steel 
to the micrometeoroid environment and measures penetration rates as a 
function of material thickness. These experiments fulfill an immediate 
need for design information. 

In the second type of experiments, which may be called scientific 
experiments, one attempts to obtain space environment information 
which, ideally, does not depend on the instrument used. In the case 
of micrometeoroids one would be interested in the spatial distribution, 
the variation of mass with abundance, and the velocity distribution; 
as well as density, velocity, and mass of individual micrometeoroids. 

So far the only operational instrument in the scientific class is the 
momentum detector which measures the product of mass and velocity of an 
impact. It is clear that if such an instrument could be coupled with 
a velocity or an energy detector, the more fundamental physical quanti- 
ties of velocity and mass could be separated. 
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In one type of velocity measuring apparatus two plastic capacitor 
foils of known spacing were used (ref. 3)- The particle must penetrate 
the foils and the time interval which elapses between the discharges of 
the two capacitors is a measure of the particle velocity. However, to 
prevent the particle's breakup and velocity reduction such foils must 
be extremely thin. This suggests that when very small particles at 
high velocities must be detected, a method of measurement is desirable 
which does not alter any parameter of the particle and which preserves 
the particle for further analysis. In this report the physical limita- 
tions of such an instrument will be investigated. To make an estimate 
of these limitations, all the available environmental information must 
be used. 

From visual observations of meteors, radar measurements, rocket 
and satellite observations, a crude estimate of the spatial distribu- 
tion of micrometeoroids has been made. For the purpose of this report 
only the data are of interest which help to determine the design param- 
eters of the instrument. The most important piece of information is 
an estimate of the mass versus flux distribution of the meteoroid popu- 
lation. Figure 1 shows various estimates and measurements (ref. 4). 

The diagram vividly points out how much there is yet to learn. Other 
important pieces of information needed are the densities and the shapes 
of the cosmic dust particles. Data of this kind has been obtained by 
high altitude rocket studies (ref. 5)* The measurements shown in fig- 
ure 1 all have been momentum measurements. The data have been reduced 
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to mass by estimating an average particle velocity. At this time no 
measurements have been made which determine both velocity and momentum 
of individual cosmic dust particles. 

The micrometeoroid velocity detector described in this report 
consists of three main types of components, the housing, the two light 
screens, and the four photo -multiplier tube detectors (fig. 2). From 
a specified field of view the instrument housing accepts micrometeoroids 
and it unavoidably accepts also an amount of undesirable light. 

By proper geometrical design and use of a special black paint, an 
attempt is made to reduce stray light falling on the face of the photo- 
multiplier tube. 

The light screens are rectangular beams of focused sunlight. 

Focusing is necessary to achieve the light intensity needed to detect 
the small size particles. 

When a cosmic dust particle traverses the first light screen, light 
is scattered and reflected from the particle and a small fraction of 
this light reaches the cathode of the photomultiplier tube. The result- 
ing current forms the start pulse of the detector. A similar output 
from the second screen forms the stop pulse, and the time elapsed between 
the two events is electronically measured. Since the distance between 
the screens is known, the time of flight between the screens determines 
the velocity of the dust particle. The addition of a momentum detector 
behind the light screens allows the determination of particle mass. 

I 
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The principle of operation of the cosmic dust velocity detector 
is very simple. Still a surprising number of problems must be solved. 
Only a few of these are treated in detail in the following report 
which may be considered a feasibility study. For completeness, some 
of the problems not treated later are enumerated here, (l) The pro- 
duction of high-speed small particles for testing of the instrument has 
not been solved satisfactorily. (2) No consideration is given in this 
report to the actual circuit design. (3) Power, weight, and tempera- 
ture control considerations in space are not treated. (4) The method 
of producing the light screens by proper focusing of the sunlight and 
the difficult problem of stabilizing the instrument to maintain this 
focus are not investigated. 


SYMBOLS 


a 

A r 

C 

d 

e 

E 

E coll 


width of the light screen 

relative front window size, (number of particles passing 
through front screen per day)/(pumber particles striking a 
one m 2 area per day) (See fig. 1 for the denominator.) 
photomultiplier tube load capacity 
distance between light screens 
electronic charge 1.602xl0~ 19 coulomb 
error rate of the velocity detector 

collection efficiency of the optical system, (number photons 
collected at one photo cathode^number photons scattered) 
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E quant 

E scat 

G 

i 

I 

Ic 

K 

L 

m 

n 

n l 

n 2 



n I 

n II 


quantum efficiency of the photomultiplier, (number electrons 
emitted from photo cathode) /(number photons collected at 
1 photo cathode) 

scattering efficiency of the particle, (number photons 
scattered) /(number photons incident) 
current gain of the photomultiplier tube 
current , ampere s 

intensity magnification of the sun's energy 
number of particles going through the circular area of 
radius r 

proportionality constant 

thickness of the light screen, through which the particle 
has to travel, cm 
mean of a probability distribution 

number of photons arriving at the photo cathodes of one 
light screen 

average number of errors per second when 1 additional pulse 
occurs in the interval t m 

average number of errors per second when 2 additional pulses 
occur in the interval t m 

n 1 

binomial coefficient, 

x! (n - x) I 

noise pulse rate from the first light screen 
noise pulse rate from the second light screen 


i 
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n q 


N 


av 


N e 


P(-) 

p (y/x) 


q. 

r 

R 


R 


d 


R c 


m 


T- 


m 


a 

V 


x 

y 

z 


number of light quanta arriving in the vicinity of the earth 
from the sun in the range from 2000 to 8000 Angstrom, 

4 . 8 5x10 16 ( quanta/ cm 2 s e c ) 

collection rate, average number of micrometeoroids being 
detected per day 

average number of electrons required from the photo cathode 
to get a high probability of an output distinguishable from 
noise, electrons 
probability density function 

probability that event y will occur given that event x 
has occurred 
charge, coulomb 

equivalent radius of a particle, cm 
resistance, ohms 

phototube distance from a point source of light, cm 
coincidence rate, (number of particles going through both 
screens per dayj/number of particles going through front 
window only) 

velocity of the particle, cm sec -1 

maximum flight time of a particle through a light screen, sec 
maximum flight time of a particle between light screens, sec 
peak voltage of a photomultiplier output pulse due to one 
primary electron, volts 

number of primary electrons due to n photons 

minimum number of output pulses to be considered an event 

standardized normal variable 
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a, (3 direction angles, defined in figure 5 

a average number of noise pulses per second 

p target size reduction factor 

a variance 

dE differential solid angle 

r pulse ■width 

$> influx rate, number of micrometeoroids above the detection 

threshold striking aim 2 area per day (See fig. 1.) 

$ 0 light flux from a zero magnitude photovisual star, 

2.43X10 -10 lumen cm -2 

Oiji total starlight flux through a unit area, lumen cm -2 

THE SENSITIVITY OF THE MICROMETEOROID VELOCITY DETECTOR 

The power of resolution of the micrometeoroid detector is defined 
as the minimum equivalent particle radius that the instrument can 
detect with some certainty. The equivalent particle radius is the 
radius of a circle which has the same area as the projected area of 
the micrometeoroid (see fig. 3)* The degree of certainty of detection 
that is desired depends on the number of particles which traverse the 
instrument. For instance, if several hundred particles above a certain 
size traverse the instrument, 80 -percent detection would be sufficient 
to infer as to the total number of particles which actually did traverse 
the instrument. The design of the transducer, which includes the light 
screen, the light collector, and the photomultiplier, offers the most 
serious limitation to the size of the particle that may be detected with 
any degree of certainty. 
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If the previously defined quantities are used, the equation for 
minimum, detectable particle size will look as follows: 


r(cm) 


1 111 1 1 1 

L n q I E sca ^ E co p]_ Eq Uan ^- 



( 1 ) 


By writing a unit equation the above equality may easily be checked. 

The individual factors of this equation will now be discussed further. 

A diagram of the detector's light screen arrangement is shown in 
figure 2. The thickness of the light screen, L, through which the par- 
ticle has to travel is limited by the size of the instrument and also by 
the requirement that stray light must be kept to a minimum, since stray 
light constitutes noise. The velocity range of the particles is esti- 
mated from the following astronomical constants. The escape velocity 
from the surface of the earth is 11 kilometers per second. This would 
be the lowest velocity one would expect cosmic dust to have in the 
vicinity of the earth. The highest velocity of micrometeoroids 
expected is 72 kilometers per second. This figure is derived from the 
following constants. The escape velocity from the sun's gravitational 
field at one astronomical unit is 42.2 kilometers per second while the 
earth velocity around the sun is 29-8 kilometers per second. Hence, 
the maximum collision velocity of a micrometeoroid bound to the solar 
system with the earth would be 72 km sec -1 . Particles with speeds 
greater than 72 km sec -1 have velocities of hyperbolic orbits. Such 
particles are not likely to be intercepted as they are probably smaller 
in number. Satellite motions with respect to the earth have been 
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neglected in the above discussion. They will increase further the 
expected spread of velocities. 

A great amount of theoretical information is available on light 
scattering of spherical particles of different sizes and different 
indexes of refraction. No theory, however, exists for light scattering 
of irregular shaped small particles . One thing that can be learned 
from the plotted scattering distributions for small spherical particles 
is that the intensity for different scattering angles will vary widely 
with wavelengths (ref. 6). This points out the fact that it would be 
undesirable to use single frequency light as that of a laser, for 
instance, even though the use of a single frequency light could provide 
a method for keeping stray light from the detector by filtering the 
incoming light to the phototube. The relatively wide band width that 
sunlight provides helps to randomize the scattering angles and thus 
makes the photons arriving at the photo-cathode more likely proportional 
to the actual projected area of the particle. 

In June 1961, the Venus Fly Trap experiment actually sampled 
micrometeoroids from an altitude of 88 to 168 kilometers above the 
earth surface (ref. 5). Basically, three types of micrometeoroids were 
found; nearly spherical particles of 2 to 3 microns in diameter, 
irregular specimens which almost have the appearance of museum meteorites 
except for their smallness, and extremely irregular pieces of fluffy 
material. The regular spheres comprised the smallest percentage of 
micrometeoroids found in this experiment. For this reason not much 
more can be done than to estimate the percentage of photons scattered 
and to assume uniform scattering. 
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As equation (l) shows , the greater the photon flux in the light 
screen, the greater the sensitivity of the velocity detector. To 
provide the intense light which is needed to get a reasonable detector 
sensitivity it is most economical to use focused sunlight directly. 

Using 10 -percent efficient solar cells and 10-percent efficient light 
sources would give an over -all efficiency of only 1 percent, while 
with a mirror collection system a much higher efficiency can be achieved. 
The difficulty with using mirrors, however, is that very precise attitude 
control on the light collecting apparatus must be maintained. This calls 
for a complex instrument; however, the use of such an instrument is not 
ruled out as satellites get bigger and more sophisticated. 

The collection efficiency of the velocity detector is a function 
of its geometry. For a very simplified model of the telescope and its 
input see figure 4. The micrometeoroid is a point source radiating 
equally in all directions, and the photocathode of the phototube is a 
distance away from the point source. Then the intensity of light 

falling onto the photocathode is proportional to l/R^ 2 . From this 
point of view it is clear that it would be desirable to have the photo- 
multiplier tube as close to the micrometeoroid as possible. However, 
there are two more considerations. Light from various sources enters 
through the same opening through which micrometeoroids enter the instru- 
ment. This light constitutes noise. Also, the closer the photocathode 
is to the light screen, the greater are the variations in telescope 
sensitivity as function of the point through which the micrometeoroid 
penetrates . Figure 4 shows a better detection uniformity when two 
phototube outputs are paralleled. 
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Three principal sources of disturbing light are present: starlight, 

moonlight, and earthlight. The telescope is assumed to be protected 
against sunlight since its collection mirror mast always be sun oriented. 
Therefore, the mirror shields the telescope from the direct rays of the 
sun. 

The data for the following calculations are taken from reference 7. 
Starlight from the whole sly is equivalent to 490 zero photo -visual 
magnitude stars (ref. 7)* If the total light flux is assumed to be 
equally distributed over the whole sphere (isotropic field), the light 
flux going through one side of a flat surface of 1 cm 2 in one direction 
will be 

$T = 490 $ o = 2. 98x10" 8 lumen cm -2 (2) 

Total number of quanta due to starlight through 1 cm 2 = 2 . 98 x 10 "® lumen 
cm" 2 4. 12x10 15 quanta sec" 1 lumen" 1 = 1.23x10 s quanta sec" 1 cm" 2 . 

Satellites are close to the earth in relation to the earth-moon 
distance. The illumination at earth's surface due to full moon is 
3x10 " 6 that of the sun (ref. 8). The sun radiates 4. 8 5x10 16 quanta cm -2 
sec" 1 , the moon 4.85xl0 ls x3xl0" s = 1.455X10 11 quanta cm" 2 sec" 1 ; that is, 
moonlight at full moon is about 1000 times as strong as starlight. 

The Bond Albedo is the ratio of total light reflected from a sphere 
to total light incident on it (ref. 7 )- The earth and moon albedo are: 
earth = 0 . 34 , moon = 0.07. From these figures it is obvious that for a 
satellite in the vicinity of the earth, the earthlight is the strongest 
contributing factor to noise due to stray light. 
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Which type of light we must guard against depends on the mission 
that is flown. On earth satellites, all three sources of light are 
important. If the design is such that only noise due to starlight is 
small enough not to cause an excessive amount of errors, one can guard 
against moonlight and earthlight by not pointing the instrument in the 
direction of the moon or the earth. For space probe missions only star- 
light and sunlight need to be considered. To prevent operation under 
unfavorable conditions the output pulses of the photomultiplier tube 
from the first light screen of the detector can be counted. When the 
counting rate goes too high, indicating a noise condition against which 
the coincidence circuits cannot discriminate, the velocity detector 
will be disabled. Also, the noise pulse count can be telemetered to 
ground at all times as an indication of confidence for the recorded 
data. 

The characteristics of the phototube detector and the quantum 
nature of light will offer the most severe limitations to the detection. 
To design an optimum detection system, the characteristics of the photo- 
tube must be clearly understood. Therefore, a later section of this 
report is devoted to that topic. At this time, let us state briefly 
that a detection system is needed which can distinguish the least 
possible number of photons falling on the photocathode as an input as 
compared to random thermal and starlight noise. These requirements 
indicate a phototube with a low thermal emission and a high quantum 
efficiency, and an electronics detection system following the phototube 
which recognizes a minimum number of primary electrons as a signal 
output . 
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GEOMETRIC FACTORS DETERMINING THE USEFULNESS OF THE 
MICROMETEOROID VELOCITY DETECTOR 

Even with coincidence circuitry, due to noise, a certain number 
of false particle counts will he made. This noise counting rate must 
he at least an order or magnitude lower than the actual particle rate. 
Hence, an estimate of the number of particles that will he encountered 
must he made to determine whether a certain size velocity detector is 
of practical use. 

Then 

N av = A r R c O (counts per day) (3) 

where the quantities A r and R c must he calculated from the geometry 
of the detector and 0 must he estimated from figure 1. 

It is intuitively clear that in an isotropic field, A r must he 
equal to the ratio of the respective areas 

^ _ area of the front screen ^ 

r - area of the reference surface 

Also, this value should not depend on the shape of the screen. To cal- 
culate the average error in velocity measurement one must also know the 
directional pattern of the screen surface as a function of the angle a. 
We will obtain this information hy first considering the circular front 
window of figure 5(a). 

The average number of particles arriving from a solid angle is 
proportional to the magnitude of the solid angle. Hence, d£, the 
differential flux from direction a (for (B = 0 to 2jt) is 

lb 
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d£ = K2jtr 2 sin a da 


( 5 ) 



which is proportional to the area of the screen as was expected. 

We must now generalize to other surfaces. Figure 5(h) shows an 
area and two different projections of the same area viewed from an 
angle a for two different angles (3 . See figure 5(a) for the defini- 
tions of the angles. Since in each case the respective differential 
areas are multiplied by a constant factor cos a, the size of the pro- 
jected area is independent of the angle 3, even though its shape varies 
widely. 

Ag = A 3 = A]_ cos a ( 9 ) 

We are now prepared to generalize. Taking the results from equa- 
tions (8) and ( 9 ) one can see that (4) holds for any shape of screen. 

The expression of equation (7) can be normalized by dividing through 
by the proportionality constant and the area of the screen. 
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d[Ic(a)] 0 . 

= 2 sin a cos a da 

Kkr2 


( 10 ) 


We rewrite equation (10) by use of trigonometrical identities and define 
it as 

p(a)da = sin 2a da (ll) 

The expression p(a) of equation (ll) can be thought of as a probability 
density function since the isotropic field is of statistical nature, 
and 

p(a)da = 1 

From equation (9) p(a) is the probability density function for any flat 
surface (see fig. 6). This indicates that on the average only 3 percent 
of all particles arrive from angles between 0 - 10° (from the front), 
while 17 .5 percent will arrive from angles between 40° - 50° • The 
relative flux intensity on a surface is redrawn in figure 7 polar 
coordinates to show the directional pattern of a flat surface. The 
result is, at first glance, somewhat surprising since one would expect 
a = 0° to be the favored direction (p is a maximum at a = 0) . How- 
ever, the differential flux is a function of d£ also, which has its 
maximum at a = jc/2. 

By adapting the method used in reference 9 , the coincidence ratio 
was calculated for light screens with a cubic geometry. A value of 
0.213 was obtained for R c . This problem actually involves quadruple 
integration, resulting in a rather lengthy method when one intends to 



16 


investigate various detector geometries. For these reasons it is useful 
to have an experimental technique to measure the coincidence ratio 
directly. 

The basic difficulty lies in simulating an isotropic field of 
particles. This was solved by an experimental analog and no mathematical 
proof is given. 

A 5 *25 -inch radius plastic hemisphere was spray painted white 
inside and outside to get good light diffusion. The opening simulating 
the front screen was placed at the center of the hemisphere, and was 
made much smaller than the radius of the hemisphere. This limits the 
variation in light intensity on any point of the screen due to illumina- 
tion from a particular spot of the sphere's surface. The remaining 
error is largest when the light comes from a direction close to the 
plane of the front screen, but this effect is minimised due to the 
smaller projected area of the screen. 

The hemisphere is inserted into a large box painted white on the 
inside (see fig. 8) where ten small light bulbs are distributed to 
illuminate the sphere with approximately equal intensity from all 
directions. 

The rest of the system is explained in figure 8 . The variac is 
used to increase the light intensity as the distance between front 
screen and back screen is increased in order to stay above the noise 
level of the tube. At points where the light intensity is increased, 
two measurements must be made. 
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The coincidence ratio is simply the ratio of photomultiplier output 
voltage measured at distance d (corrected hy any light increase factors) 
to the voltage measured when d = 0. 

The results are plotted in figure 9* The result for d = a which 
was computed earlier is very close to the measured value. Notice the 
rapid decrease in coincidence ratios for small separations between 
screens. As d increases further, the front face appears as a point 
source and further decrease of light intensity versus distance is 
proportional to l/d 2 . 

Several geometric arrangements of the telescope can be investigated, 
including a shield of some sort which will be needed to keep out star- 
light and reduce the noise background. One type of shield is shown in 
figure 10(b). Less particles than shown in figure 10(a) will now go 
through the front screen, but the number of particles going through the 
back screen have not been reduced since the shield does not cut off any 
coincident rays. Hence, the collection rate of coincident particles 
can be computed as if the shield were absent. The computation would not 
indicate, however, the reduction of particle flux through the front 
window. If this is of interest, special measurements must be made. 

Figure 10(c) shows a straight shield. Obviously this type of 
shield reduces the number of coincident particles captured when com- 
pared to figure 10(a). See table I for a comparison of the perform- 
ances of the different detector geometries. 

The above method can be extended to find the acceptance pattern 
of any velocity detector configuration by the simple method of masking 
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different parts of the hemisphere and recording the percentage of light 
received from the unmasked portion of the hemisphere. The acceptance 
pattern must he known for two purposes, to determine the errors in 
velocity measurement because the particle does not traverse the shortest 
flight path, and to determine the direction of the measured micrometeoroid 
flux. 


PHOTOMULTIPLIER TUBE CONSIDERATIONS 

In the micrometeoroid velocity detector we are looking at small 
particles passing a 1-centimeter illuminated distance at speeds between 
7 and 72 kilometers per second. Hence, the duration of events to be 
detected ranges from 0.14 to 1.4 microseconds. From equation (3) and 
figure 1 using Whipple's data, it is estimated that on the average about 
four particles of a 2-micron radius would be detected per day. The same 
instrument would detect on the average a 100 -micron particle every 
1000 hours. This indicates that to get a statistical sample the photo- 
multiplier tube must be used in an optimum fashion to detect a minimum 
number of photons arriving at the photocathode. Since the large gain 
of the tube makes it possible to distinguish 1 photoelectron, it becomes 
necessary to investigate the noise behavior of the tube very carefully 
(ref. 10). 

To 'understand noise behavior of a tube, the gain of the tube must 
be known. Measurements were made with an EMTQSOP'B* tube which has the 
following nominal characteristics: For an over -all voltage of 

*Electrical and Musical Instrument Co., Los Angeles, Calif. 
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1500 volts across a uniform dynode chain it has a sensitivity of 
2000 amperes per lumen. From this information the current amplifica- 
tion can he estimated: 


2000 ■ 0.624x10 10 electron/sec amp 

G = 1HSS2 ' = 2-75X10 7 

4.53X10 15 g2- o.i electro11 
sec lumen quantum 


Another way of arriving at an approximate value of current amplifica- 
tion is found by estimating the dynode multiplication factors. These 
are equal to 6 for the amplification from the photocathode to the first 
dynode and equal to 4 between all other dynodes. Hence, for the 
13 -dynode tube we calculate: 

G = 6x4 12 * 10 s 


Since the approximate gain of the tube was known, the actual tube gain 
was measured more easily (see fig. 11). 

Figure 12 shows a circuit to measure the dark current pulse rate. 

The usual calculation for the mean pulse height is 

''ideal = c = f (13) 

For figure 12 

$ = 3-3xlo 7 xi.6xio- 19 _ 0 -1(6t 

11. 5x10" 12 

This assumes a current impulse is delivered, so that current through 
the resistor during the charging of the capacitor is negligible. When 
this is not the case, the pulse height will be lower than that calculated 
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above. Figure 13 shows a Gaussian current pulse. For this actual 
pulse shape it is difficult to calculate the output voltage. There- 
fore, let this pulse he roughly approximated hy a rectangular pulse 
having the same total charge, and a width equal to the half -amplitude 
time spread (i = q/r) . Then the actual peak output voltage for a 
pulse is reduced hy a factor 


A 


^actual 

^ideal 




(14) 


due to the presence of resistance. In figure 14 the relationship of 
equation (l4) has heen plotted. It shows that for a given load capaci- 
tance the pulse height increases with increased load resistance. 

For the EMI9502B with a load as shown in figure 15, 


RC = 1,1. 5X10-3 = q>64 

T l8xl0" 9 

^actual = V id eal( RC / T )( 1 - e- T/RC ) = 0 . 46x0 . 64x0 . 98 =0.293 volts 
is the average output pulse height due to 1 photoelectron. This value 
has heen used as a scale in drawing the abscissa of figure 15 - 

Since thermal electron emission is not a function of supply voltage, 
the quantity eG is useful as the scale for the abscissa because fig- 
ure 15 remains the same even when the tube gain is altered by changing 
the supply voltage. To obtain an estimate of the proportion of pulses 
originating at the photocathode, the photocathode was reverse biased 
by 100 volts. The pulse rate for small pulses (0.33 eG) decreased by 


21 



a factor of 10, and the pulse rate for larger pulses (0.66 eG) decreased 
hy a factor of 150. Hence, nearly all large pulses originate at the 
photocathode, -which is not a surprising result. 

When coincidence methods are to he used to reduce counting errors, 
it is important to know the time distribution of the pulses (ref. 11). 
Figure 12 shows the equipment required to measure this distribution. 

CRO II was free running at a known speed and the scope face was masked 
except for 1 centimeter of height. For each pulse above a minimum 
height CRO I swept once and applied a pulse to the slow sweeping CRO II. 
A photograph of the face of CRO II is shown in figure l6. By consider- 
ing the last event of the first row identical with the first event of 
the next row (even though they are not) we have a relatively long 
sample of pulse versus time distribution. 

One can now test for the assumption that the pulses form a Poisson 
time distribution. This is a reasonable assumption to be tested since 
such a distribution will result provided the following conditions hold: 

1. The probability that n electrons are emitted within the 
interval (t 0 , t 0 + t) is independent of the starting time t 0 . 

2. For an infinitesimal interval of length, At, the probability 
that 1 electron occurs within At is proportional to At, and the prob- 
ability that more than 1 electron occurs within At is negligible. 

3. The time of emission of individual electrons within intervals 
that do not overlap is statistically independent. 
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Then the probability of exactly n electrons being emitted in the 
interval t is : 

Pn = P(n electrons in the interval t) = ^ - a ^ — e" a ^ (15) 

A n! 

where at is the average number of pulses in the interval t. 

One can now count the number of pulses in figure 1 6 to determine 
the average number of pulses per second. By dividing figure 1 6 into 
intervals of length t = l/a, one can count the number of pulses in 
each interval and make a table of their distribution. See table II 
and note the agreement between theoretical and measured values. 

Environmental temperature drastically affects the noise output of 
the tube. With decrease in temperature the photomultiplier tube gain 
remains virtually unchanged, while the noise pulse rate is divided in 
half for every 10° C temperature decrease. Therefore, pulse height 
versus pulse rate was measured again when the tube was cooled to -50° C 
and it was found that the Poisson law still applied. 

What does this mean in terms of the micrometeorite detector? If 
the maximum temperature at which the tube would work on the satellite 
is known, the pulse height distribution of the tube for this tempera- 
ture can be measured, and the one parameter noise statistics of the 
tube would be completely described. 

Unfortunately, when the photomultiplier tube is integrated as part 
of the micrometeoroid detector, light from the stars will arrive at the 
photocathode and thus constitute noise. At first glance, it seems that 
this light should cause a steady-state dc value. However, for very 
weak light inputs, this is not so. 
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To determine experimentally what happens , the tube was placed in 
a container which was designed to let in a minimum amount of starlight 
and still allow the system to perform as a micrometeoroid velocity 
detector (see fig. 2). The instrument and its associated measuring 
equipment were setup on Mount Hamilton, California, during the new- 
moon period in August 1962. The experimental set-up is shown in figure 17 • 

The results of the experiment are illustrated in figures 18, 19, 
and 20. As figure 18 shows, the pulse shape due to the starlight noise 
is identical with that of thermal noise alone. However the frequency 
of the noise pulses increased by 1-1/2 decade (fig. 19). With such an 
increase, extreme cooling of the phototube would give rise to little 
improvement in the over-all performance of the detector. It is hoped, 
however, that eventually a considerable reduction in noise due to star- 
light could be achieved by more careful experimental design of the 
,f b lack -box" structure. 

The statistical characteristics of the noise pulses were tested 
by the same method as described earlier and were found also to be 
Poisson (see fig. 20). This follows because the addition of two Poisson 
processes will result in another Poisson process with a new parameter 
a being the sum of the two individual parameters. 

THE SYSTEM 

In a detector two types of errors can be made. The velocity of a 
micrometeoroid going through both light screens would not be measured 
as a result of the failure of an output from at least one of the light 
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screens (this is called a type I error); or, as a result of the noise, 
hoth screens could give an output in the proper sequence when actually 
no event has taken place (this is called a type II error) . Designing 
the experiment to minimize one type of error will necessarily increase 
the prohahility of the other type of error; therefore, an acceptable 
compromise must be found. The object is to count a large proportion 
of all particles entering the detector and to get an average error 
count of only a few percent of the actual count rate. 

Errors of type I are estimated in the following manner. Let n 
be the number of photons arriving at the photocathode, and p be the 
quantum efficiency of the phototube. Then the probability that 
exactly x photoelectrons are emitted from the photocathode is a 
binomial distribution (ref. 12) 

P(x) = Q P x (i - p) n ' x (16) 

with mean and variance 

m = np (17) 

c = Jnp(l - p) (18) 

For small p (10 percent) and large n the binomial probability dis- 
tribution can be approximated by the normal distribution with the same 
mean and variance. The evaluation using the normal distribution approxi- 
mation can be shown to be (ref. 13) 
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where the values of the standardized normal variable 


P r (Z < Z«) = — L f e" z2/2 dZ (20) 

n/Sj t '-'-oo 

can he determined from standardized normal tables. 

It was shown earlier that only one -tenth of the pulses which are 
smaller than one -third of the average output pulse height originates 
at the photocathode. This indicates that one can reduce type I errors 
considerably by biasing the pulse discrimination circuits that follow 
the phototube output. To suppress virtually all pulses not originating 
at the photocathode, the discrimination circuits may be biased to dis- 
tinguish pulses greater than(l/2)eG. The probability that the number 
of output pulses N is greater than or equal to y given that x 
primary electrons are present is 
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Since the processes expressed in equations (19) and (21) are independent, 
then (ref- 10) 

P(x,N > y/x) = P(x)P(N > y/x) (22) 

Thus the probability of an output when n photons arrive at the photo - 
cathode is 

n 

Output = X p(x)p(n > y/x) (23) 

n=y 

This sum can easily be evaluated for a different number of incoming 
photons n and for a different minimum number of output pulses y which 
are considered an event. A sample calculation is shown in table III. 
Before discussing the significance of these figures let us estimate 
type II errors. 

To reduce type II errors, coincidence methods are usually used when 
outputs from separate detectors axe involved (ref. 10). However, delayed 
coincidence can also be applied when successive outputs from a single 
detector concerns us. 

One can assume the following circuit. It gives an output pulse 
when, after an initial pulse, a minimum number of additional pulses occur 
within a given time. This time is chosen as t m the maximum flight time 
of a particle in the light screen. 

As previously shown the noise pulses occur according to the Poisson 
distribution 
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P n (t m ) = P n (n noise pulses in the interval 0 < t < t m ) 

_ (cxtni) e -at m (2b) 

nl 

•where a is the average numher of photomultiplier noise pulses per 
second. After each pulse the circuit counts for t m seconds if the 
specified minimum number of pulses occur, for example, one extra pulse: 

P n (n > 1 in 0 < t < t m ) =1 - P(0) 

= i _ e “ atm 

* 1 - (l - at m ) for at « 0.1 
= at m (25) 

or two extra pulses 

P n (n > 2 in 0 < t < t m ) = 1 - P(0) - P(l) 

-atm fc^m. 0 (atm) 

" L 01 1: 

~ 1 - (l - at m ) ( 1 + atm) 

= (at m ) 2 (26) 

There are an average of a counting intervals per second; hence, 
the probability that in these a counting intervals an event is indicated 
falsely is 

aP(n > 1) = a(atm) = cc. 2 t m = n x (27) 

aP(n > 2) = a(at m ) 2 = a 3 t m 2 = n 2 (28) 

where n ± and n 2 are the average number of errors per second; for 
example, a = 1000 pps, t m = 2[isec 
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ri! = ( 1000) 2 2x10 " 6 = 2 errors per sec 
n 2 = ( 1000) 3 4x10" 12 = 4x10" 3 error per sec 


In the above circuit it was assumed that the individual noise pulses 
are short when compared to t^. A O.lpsec pulse can reasonably be 
achieved, and indeed t m = 2|j.sec » O.lpsec. 

Clearly a single coincidence circuit with Pa 1 ( r i > 1) or even 
P a 2 (n > 2) cannot sufficiently reduce the effect of noise, when rare 
events must be detected (10-200 events per day). To reduce type II 
errors further, another coincidence arrangement which involves both 
screens can be employed. Again, the front screen mast first indicate 
an event; therefore, the equation for noise improvement is identical 
with that previously developed. 

Figure 21 shows that a 20|-isec pulse is applied to the coincidence 
gate every time the first screen gives an output. This allows a maxi- 
mum flight time between screens of 20p.sec and slower particles will not 
be detected. From the second screen a 2|_isec pulse is sent and checked 
for coincidence. In terms of the noise in the detectors this results 
in the following; 


E 


errors 

seconds 


- T m n T n- 


I U II 


(29) 


where nj and njj are the noise pulse rates from the first and second 
screens . Since T m = lOt^ (the screen distance is 10 times the screen 
thickness), the following holds 


E = 10t m njnjj 


(30) 



Two cases have "been considered which resulted in n x and n 2 errors per 
second (eqs. (27) and (28)). 

If noise pulse rates are assumed to he equal for both velocity- 
screens, then the first case (eq,. (30)) will become 

E = 10t m 3 a 4 (sec _1 ) (31) 

which gives 

a = 4/ — "~~3 = 1.06xl0 4 ^fE (32) 

kJ lOtm 

For the second case 

E = 10t m s a 6 (sec _1 ) (33) 

which gives 

a - J35 - 3 - 82xio4?fs m 

The permissible noise pulse rates for one error per week and one error 

per month are shown in table V. It clearly shows the superiority of 
the second method. 

A block diagram can now be drawn and the choice of logic circuitry 
justified (see fig. 22). For performance calculations the following 
assumptions are made: 

1. If the geometry in figure 22 is used, starlight will contribute 
to the tube noise to give a maximum number of pulses of 2000 per second 
per tube. 
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2. Tubes are operated at -50° C. At this temperature the 2-inch 
diameter photocathode will emit about 6 thermal electrons per second, 
which is a negligible amount compared to starlight noise. 

3. Output pulse length from one photoelectron (noise or signal) 
is O.lpsec for both tubes in parallel monitoring the same light screen. 

4. Sunlight in the light screen is intensified by a factor of 100. 

5 . Particles scatter photons equally over 4it steradians for the 
calculation of collection efficiency. 

6. Scattering efficiency of the particles is 10 percent. 

7 . Average particle velocity is 15 km/sec (from Whipple) . 

8. Phototube quantum efficiency is 10 percent. 

9 . Particle influx rate follows the curve given in reference 4 

and figure 1. (This curve gives a higher influx rate for small particles 
when compared to Whipple's 1957 data.) 

10. Three output pulses within 2psec from one phototube pair are 
considered an event. 

Calculated Performance Data 

1. The probability of detecting the smallest size particle considered 
is 78 percent (see table IV). 

2. The equivalent radius of the smallest particle seen with the 
above reliability is from equation (l). 


r 


1 1 _ 1 _ 1 1 1 1 

it L n q I E sca t E co q]_ Eq uan ^ 




1.5X10 S 


1 1 _1 1 _ 1 1 _ 

4 4.85X10 10 10 2 10' 1 1.8X10" 2 10" 1 



= O .73 micron 
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Above this size the probability of detection rapidly approaches 1. 
Below this size the probability of detection rapidly approaches zero. 
This can be verified by the method of calculating type I errors. 

3. From equation (3) and figure 1 (McCracken ) , the average number 
of 2. 5 -micron particles counted per day in the vicinity of the earth 
is 


N a v _ 


100 


m 2 0.2 1 


particle 
m 2 - sec 


^ particle _ particles 
500 sec day 


(Whipple’s estimate would result in a figure of N av = 4.4 particles 
per day . ) 

4. From equation (33) the average number of errors due to noise 
is 

E = 10 ( 2x10 ” 6 ) 5 ( 2xl0 3 ) 6 

= 2x10-8 errors 
second 

= 1 error in 0.63 year 

This error rate is negligible, but one must consider the fast increase 
of errors as noise pulses increase (see table V) and be aware that the 
margin is not as great as it seems . Further improvement can be achieved 
by integrating the momentum detector into the coincidence scheme. 
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CONCLUDING REMARKS 


It appears that special techniques would he required in order 
to achieve detection of micron-sized particles by light screens. 

The suggested method involves the sensing of pulses due to individual 
electrons originating at the photocathode of a photomultiplier tube. 
The analysis of the two basic types of errors, failure to detect a 
particle, and erroneously interpreting noise as a particle shows that 
a coincidence type of detection scheme can make the latter type of 
error reasonably small. Even after going to these lengths to avoid 
errors due to noise, an extremely intense light (of the order of 
100 times the intensity of direct sunlight) would be required. This 
need for high intensity light screens may well be the deciding factor 
as to the practicability of the device for use in a space vehicle. 

The suggested method of detection and noise discrimination in 
which the pulses resulting from individual electrons originating from 
the photocathode are scrutinized is believed to be new and may have 
wider application. For example, it might be applied to scintillation 
counters for the detection of extremely low-level rare nuclear events. 
The analog method used to investigate the directivity and coincidence 
ratio of the velocity detector should also be useful for other types 
of detectors. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 9 > 19^3 
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TABLE I . - VELOCITY DETECTOR PERFORMANCE COMPARISON WHEN A SHIELD IS 
USED ASSUMING THAT ON THE AVERAGE 17 PARTICLES WILL ENTER A 
10- BY 10 -CENTIMETER AREA PER DAY 


Type of screen 

2 -micron particles 
entering front 
screen/day 

2 -micron particles 
entering "both 
screens /day 

1 column 3 

c column 2 

Fig. 10(a) 
(no shield) 

17.0 

3.76 

0.22 

Fig. 10(b) 

d 2 = dj** 

9-4* 

3.78 

.4 

Fig. 10(c) 
d 3 = 5 cm 

8.5 

1.88 

.22 

Fig. 10(c) 
d 3 = 10 cm 

3.76 

1.1 

• 29 


*Measurement by testing a model shield 
**a = 10 cm, d = 10 cm 


Note: If Rc 1 is used in calculating the number of micrometeorites 

counted per day on the average, the flux rate cp must be reduced 
to cp 1 accounting for the effect of the shield. That this must 
be done is clear because the detectors of figures 10(a) and 10(b) 
have exactly the same count rate. 


TABLE II . - COMPARISON OF THE CALCULATED POISSON 
AND THE MEASURED DISTRIBUTION 


h 

p - JL. p ~ i 
r n _ r e 
nl 

Measured 

0 

0.368 

0.38 

1 

.368 

•34 

2 

.184 

.21 

3 

.0613 

.06 

4 

• 0153 

.01 

>4 

.0034 

0 

I 

1.000 

1.00 


Note: The agreement between the calculated 

Poisson distribution and the measured 
distribution is remarkably close. 
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TABLE III . - PROBABILITY CALCULATION FOR DETECTING A PARTICLE 


(for n = 70 * y = 3 ) 


X 

P(x) 

P(N > y/x) 

P(x)P(N > y/x) 

0-2 

0.055 

0 

0 

3 

.045 

•125 

.0056 

4 

.068 

•313 

.0214 

5 

.115 

.500 

.0575 

6 

.147 

.656 

.0962 

7 

.158 

•773 

.1220 

8 

.147 

.855 

.1260 

9 

.115 

.910 

.1050 

10 

.068 

• 945 

.0640 

11 

.045 

.961 

.0430 

12 

.022 

.980 

.0215 

13 

.013 

.989 

.02128 

14 

.001 

• 994 

.0010 

> 14 

Negligible 


I" • 6 76 


The square of this sum. (O.676) 2 is the probability 
that both screens will "see" the particle. To show 
the trend, three more points were calculated and are 
entered in table IV. 


TABLE IV.- PROBABILITY OF DETECTING A PARTICLE FOR GIVEN n AND y 


Number of photons 



arriving at the 

n = 70 

n = 100 

photocathode 



y = 2 

( 0 . 84) 2 = O.7O 

(0.92) 2 = 0.85 

y = 3 

(0.68) 2 = 0.46 

(0.86) = O.74 
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Figure 3 *- Definition of the equivalent particle radius 



Light screen 
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Figure 4. - Detector sensitivity as a function of the particle position in the light screen. 
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Figure 15 •- Number of pulses greater than a given pulse height versus 

minimum pulse height. 





Figure 1 6 .- Pulse time distribution. 



A-30513-3 

Figure 17- - Test setup at Mt. Hamilton. 
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p, | H* S A-30513-1. 1 

u * 1 cm 

Figure 18.- Pulse shapes of the photomultiplier tube output due to 

starlight and thermal noise. 
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Starlight 



1 J*l 

— Light duct 

Cathode size 1 cm 


© — Measured with both ports open 
and box in ditferent position 



N s 

nr ^ 


\\ 

n\ v 


•Thermal noise and starlight 
from light duct. 

-Thermal noise 


-Thermal noise and starlight 
from the front. 

-Thermal noise 


Data for curves I and II taken 
at approximately the same 
time ; same for curves Handle, 


Minimum Dulse height (volts) 

Figure 19.- Phototube noise due to starlight and thermal noise using 

EMI9502B. 






Pulse vs. time « 330 pps , t= 0 . 2 ms/cm 


A-30513-2. 1 


Figure 20.- Pulse time distribution of photomultiplier tube output due 

to starlight and thermal noise. 



Figure 21.- Coincidence gate for two outputs. 
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Figure 22.- Block diagram of the velocity detector. 
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